QCD analysis of the diffractive structure function F 2 
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The proton diffractive structure function measured in the HI and ZEUS experiments at 

HERA is analyzed in terms of both Regge phenomenology and perturbative QCD evolution. A new 
method determines the values of the Regge intercepts in "hard" diffraction, confirming a higher 
value of the Pomeron intercept than for soft physics. The data are well described by a QCD analysis 
in which point-like parton distributions, evolving according to the DGLAP equations, are assigned 
to the leading and sub-leading Regge exchanges. The gluon distributions are found to be quite 
different for HI and ZEUS. A global fit analysis, where a higher twist component is taken from 
models, allows us to use data in the whole available range in diffractive mass and gives a stable 
answer for the leading twist contribution. We give sets of quark and gluon parton distributions 
for the Pomeron, and predictions for the charm and the longitudinal proton diffractive structure 
function from the QCD fit. An extrapolation to the Tevatron range is compared with CDF data 
on single diffraction. Conclusions on factorization breaking depend critically whether Hf (strong 
violation) or ZEUS (compatibility at low j3) fits are taken into account. 



I. INTRODUCTION 



It is now experimentally well established at HERA [1,2] that a substantial fraction of ep events is contributable to 
diffraction, i.e. color singlet exchange, initiated by a highly virtual photon. Starting with the pioneering theoretical 
work of Ref . [3] , the idea of a point-like structure of the Pomeron exchange opens the way to the determination of its 
parton (quark and gluon) distributions, where the Pomeron point-like structure can be treated in a similar way as (and 
compared to) the proton one. Indeed, leading twist contributions to the proton diffractive structure functions can 
be defined by factorization properties [4] in much the same way as for the full proton structure functions themselves. 
As such, they should obey DGLAP evolution equations [5] , and thus allow for perturbative predictions of their Q 2 
evolution. 

Indeed, such parton distributions are very useful to investigate the difficult and longstanding problem of the nature 
of the Pomeron Regge singularity. On a phenomenological ground, they are the basis of MC simulations like RAPGAP 
[6] and they give a comparison basis with "hard" diffraction at Tevatron, where the factorization properties are not 
expected to be valid. Moreover, the study of diffractive parton distributions is also a challenge for the discussion 
of different approaches and models where other than leading twist contributions can be present in hard diffraction. 
Indeed, there exists strong presumptive evidence that higher twist effects may be quite important in diffractive 
processes contrary to non-diffractive ones which do not require (at least at not too small Q 2 ) such contributions. In 
fact, there are models which incorporate non-negligible contributions from higher twist components, especially for 
relatively small masses of the diffractive system. One of our goals is to take into account this peculiarity of diffractive 
processes. 

The purpose of this paper is to derive parton distributions of the Pomeron from QCD fits of diffractive DIS cross- 
sections determined at HERA [1,2]. By comparison with previous determinations contained in the experimental 
papers, we introduce some new techniques of parametrizations and discuss various points which have appeared in the 
domain: 

-the precise determination of the effective Pomeron and Reggeon trajectories in Regge fits 

-the determination of the parton distributions at Next-to-Leading Order (NLO) in the Pomeron comparing two 
methods, either using cuts in order to damp the expected higher twist contributions or by a global fit including a 
higher twist contribution taken from models 

-the differences between parton distributions derived from ZEUS and HI data 



*Service de Physique des Particules, CE-Saclay, F-91191 Gif-sur-Yvette Cedex, France; Brookhaven National Laboratory, 
Upton, New York, 11973; University of Texas, Arlington, Texas, 76019 
^Service de Physique des Particules, CE-Saclay, F-91191 Gif-sur-Yvette Cedex, France 
* II. Institut fur Theoretische Physik, Universitat Hamburg, Luruper Chaussee 149, D-22761 Hamburg 
^University of Lund, Soelvegatan 14, 223 62 Lund, Sweden 

"Service de Physique Theorique, CE-Saclay, F-91191 Gif-sur-Yvette Cedex, France 



1 



-the need to evaluate precisely the gluon content of the Pomeron and its error since it has strong impact on the 
discussion of diffractive charm and longitudinal diffractive structure functions 

-the extrapolation of the diffractive structure function to single diffraction in pp collisions at Tevatron could give a 
precise evaluation of the factorization breaking. 

Our paper is structured in the following way: In section II, we show Regge fits and the corresponding determination 
of intercepts from HI (II- A) and ZEUS (II-B) data using a new method described in the appendix Al. In section 
III we use the same data but with additional cuts including a high (3 cut to determine the parton distributions of 
the Pomeron in a NLO DGLAP evolution framework. The possibility of two different fits for HI (III-A) and a 
remaining difference between HI and ZEUS (III-B) are thoroughly examined. In section IV, we perform a new global 
fit analysis by removing the high (3 cut, and using models to provide the higher twist contribution. In IV-A we check 
the stability of the determination of the leading twist parton distributions of the Pomeron. In IV-B we discuss the 
obtained gluon density in the Pomeron (for HI and ZEUS fits) and examine the corresponding predictions for the 
diffractive charm (IV-C) and longitudinal (IV-D) components. The final selection of parameters for HI and ZEUS 
fits is given in appendix A2. In section V we discuss the important consequence of these determinations for single 
diffraction at the Tevatron, as measured by the CDF collaboration [7]. We see that the conclusions on factorization 
breaking are markedly different for the HI and ZEUS determinations. Conclusions and outlook are presented in the 
final section VI, as a discussion of the introductory points. 

II. REGGE FITS AND DETERMINATION OF INTERCEPTS 
A. Determination of the Pomeron and Reggeon intercepts from HI data 

The diffractive structure function F^ 3 \ measured from DIS events with large rapidity gaps, can be investigated 
in the framework of Regge phenomenology and expressed as a sum of two factorized contributions corresponding to 
a Pomeron and secondary Reggeon trajectories. 

if (3) (Q 2 ,/W) = ./p/p^p^Q 2 ,/?) +/r/p(zp)^ 2 R (Q 2 ,/3) • (!) 

In this parameterisation, can be interpreted as the Pomeron structure function and F^- as an effective Reggeon 
structure function, with the restriction that it takes into account various secondary Regge contributions which can 
hardly be separated. The Pomeron and Reggeon fluxes are assumed to follow a Regge behaviour with linear trajectories 
ap,p(t) = ap.p(O) + ctp R i, such that 

fp/ P ,n/ P (xp) = / 2a P , E (t)-l d * ( 2 ) 

Jtcut Xp 

where \t m in\ is the minimum kinematically allowed value of \t\ and t cut = — I GeV 2 is the limit of the measurement. 

In the first step we determine the value of the Pomeron and Reggeon intercepts (q;p j r(0)) from phenomenological 
fits of HI data. In equation (I), the values of F 2 ' are treated as free parameters at each /3 and Q 2 point. The 
values of q;p ! r(0) are free parameters while -Bp,p, and taken from hadron-hadron data (a P = 0.26 GeV , 

ctp = 0.90 GcW 2 , Bp = 4.6 GeW 2 , £>r = 2.0 GcV~ 2 ), since the data are not precise enough to determine the 
t-slopes [1]. In order to avoid the region that may be most strongly affected by a non-zero value of the ratio R of the 
longitudinal to the transverse cross-sections, only data with y < 0.45 are included in the fit. 

Note that the fit with a single trajectory does not give a good description of the data (% 2 /dof = 65.5/19 with 
statistical errors only). This confirms the observations [1,8] that secondary trajectories in addition to the Pomeron 
are required to describe diffractive ep data. A much better fit is thus obtained when both a leading (P) and a 
sub-leading (1R) trajectory are considered. 

Since the sub-leading exchange is poorly constrained by the data, a second fit is performed to the HI data in 
which the values of F^(Q 2 , (3) are taken from a parameterisation of the pion structure function [9], with a single free 
normalization. We use the pion structure function of Ref. [9] to be able to extrapolate 1 our fits to lower Q 2 . In 
principle, an interference term could be introduced in equation (I), but this does not modify the \ 2 of the fit. 



1 Note that the Duke Owens parametrisation [10] used in Ref. [1] is not valid at low Q 2 , as it has a constant value for a Q 2 
less than 4 GeV 2 . 
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The resulting value of aip(O) = 1.20 ± 0.09 is in agreement with that obtained in Ref. [1] and is significantly larger 
than values extracted from soft hadronic data (ap ~ 1.08). Also, we find ap(0) = 0.62 ± 0.03 2 . The use of the pion 
structure function in our treatment is questionable as the Reggeon intercept is close to the rho value and not to the 
pion one. However, the HI data where a proton or a neutron is tagged in the final state can be described nicely by 
pion and Reggeon exchanges where the Reggeon structure function is also taken to be the pion one [12]. We thus do 
not expect a significant bias due to this assumption. 

A well-known technical problem for the determination of the Regge fit at each j3 and Q 2 point is the correct 
identification of the actually free parameters. We use a new method described in appendix Al and perform a direct 
fit of ap(0) and ap(0). Indeed, this method is designed to take into account the decorrelations between the parameters 
and to give a more reliable minimization procedure. Without using this method, many local minima are present in 
the x 2 which might give an incorrect value of the exponents. 

This method leads to better and faster convergence properties for the Regge fits procedure when applied to HI data 
since the number of parameters to be fitted is much smaller. We find the following values : ap(0) = 1.20 ± 0.02 and 
Q!r(0) = 0.62 ± 0.02 (statistical errors only), which are close to those determined previously. We keep these values 
fixed in the next sections. 



B. Determination of the Pomeron and Reggeon intercepts from ZEUS data 

We performed the same study using ZEUS data and we find a Pomeron intercept ap = 1.13 ± 0.04, compatible 
with what has been obtained by the ZEUS collaboration [2]. This value is lower than the HI value 3 extracted in 
the previous section, however, the selection of the diffractive sample in ZEUS is different from HI [2]. We will come 
back to this last point in the following. The Reggeon contribution is also found to be compatible with zero since the 
method of measurement (the so-called Mx method) is not sensitive to the Reggeon component [2] . 

III. EXTRACTION OF PARTON DISTRIBUTIONS IN THE POMERON 
A. Parton distributions from HI data 

As we have already indicated in Equation (1), the diffractive structure function F 2 D ^ can be investigated in 
the framework of Regge phenomenology. Moreover, it has been suggested that the Q 2 evolution of these structure 
functions may be understood in terms of parton dynamics [3], i.e. as coming from leading twist, perturbative QCD 
contributions where parton densities are evolved according to DGLAP equations [5]. 

This idea has been tested in Ref. [1] with an iterative method using the CTEQ evolution code [13] over infinitesimal 
steps in x space to solve the NLO DGLAP equations [5]. In the following we redo this analysis of HI data with a new 
polynomial method to solve DGLAP equations at NLO [14]. 

We assign parton distribution functions to the Pomeron and to the Reggeon. A simple prescription is adopted in 
which the parton distributions of both the Pomeron and the Reggeon are parameterised in terms of non-perturbative 
input distributions at some low scale Qq = 3 GcV 2 . 

As we said in the last chapter, the pion structure function [9] is assumed for the sub-leading Reggeon trajectory 
with a free global normalization to be determined by the data 4 . 

For the Pomeron, a quark flavour singlet distribution (zS q (z, Q 2 ) = u + u + d + d + s + s) and a gluon distribution 

(zG(z,Q 2 )) are parameterized in terms of coefficients and cj G ' at Q 2 — 3 GeV 2 as it was done in Ref. [1] such 
that 



introducing interference terms reduces the value of the Reggeon exponent to 0.5, compatible with regge phenomenology. The 
fact we do not put any interference term leads to an effective value of the Reggeon exponent of 0.62. On the other hand, the 
value of the Pomeron exponent remains constant with or without interference. 

3 This last value is more compatible with soft physics if one considers a global fit to all data from np, Kp and pp including the 
1800 GcV pp cross-section [11]. 

4 We checked that changing the pion structure function by some amount (20%) does not change significantly the parton 
distributions. 
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zS(z,Q 2 = Q 2 ) 
zG{z,Q 2 = Ql) 



3 = 1 



3 = 1 



(3) 
(4) 



where z = x i / ]P is the fractional momentum of the Pomeron carried by the struck parton (in the following, we use 
indifferently z or j3), Pj(C) is the j th member in a set of Chebyshev polynomials, which are chosen such that Pi = 1, 
P 2 = ( and P j+1 (0 = 2CPKC) - Pj-iiQ. 

A sum of n = 3 orthonormal polynomials is used so that the input distributions are free to adopt a large range of 
forms for a given number of parameters. Any bias towards a particular solution due to the choice of the functional 
form of the input distribution is therefore minimized. 

The trajectory intercepts are fixed to ap = 1.20 and or = 0.62 5 (as explained in the last section). Only data 
points of HI with Q 2 > 3 GeV 2 , (3 < 0.65, M x > 2 GeV and y < 0.45 are included in the fit in order to avoid large 
higher twist effects and the region that may be most strongly affected by a non-zero value of R. We will discuss this 
point in the following when we take into account the higher twist contribution to the diffractive structure functions. 
The fits include 161 data points for 7 parameters (3 for the sea quark density, 3 for the gluon density and 1 for the 
normalization of the Reggeon contribution). Thus, we have 154 degrees of freedom. 

To get solvable evolution equations, the parton distribution functions must approach zero as z — ► 1. This is 
achieved by introducing the exponential term with a positive value of the parameter a. Unless otherwise indicated, 
in the following fits a is set to 0.01 such that this term only influences the parametrisation in the region z > 0.9. 
This term is only present to ensure the convergence and plays a role in a domain where we do not include data. (For 
instance, at [3 =0.65, its value is equal to 0.97.) 

The functions zS and zG are evolved to higher Q 2 using the next-to-leading order DGLAP evolution equations. 
The contribution to F^(f3,Q 2 ) from charm quarks is calculated in the fixed flavour scheme using the photon-gluon 
fusion prescription given in [15]. The contribution from heavier quarks is neglected. 

No momentum sum rule is imposed because of the theoretical uncertainty in specifying the normalization of the 
Pomeron or Reggeon fluxes and because it is not clear that such a sum rule is appropriate for the parton distributions 
of a virtual exchange. After fitting, however, we observed that the sum rules were fulfilled within 10% despite being 
not required in the fits. 

The resulting parton densities of the Pomeron are presented in Fig. 1. One possible fit quoted here as fit 1 shows 
a large gluonic content. The quark contribution is much smaller compared to the gluon one. An other much less 
favoured fit, quoted here as fit 2, has a peaked gluon distribution at high (3 6 . It should be noticed that only the first 
solution is found if one imposes in the fit that the parton distributions have no zero over the full (3 and Q 2 range, 
whereas the second fit appears when this condition is removed. 

The appearance of a second solution in our fit to the HI data is very reminiscent of the results presented in Ref. 
[16]. In fact, the general characteristics of the two solutions obtained in the present paper and of those found in 
Ref. [16] are quite similar: the difference in the gluon content at large z (or (3), and the instability of the high-gluon 
solution with respect to cuts in the data used for the fit. Nevertheless, one has to be careful in comparing the results 



? D(3) 



the 



of Ref. [16] with the present ones: whereas the former ones represent fits to the diffractive cross section F^ 
present analysis describes results for diffractive parton densities. 

The result of the fit (fit 1) is presented in Fig. 2 together with the experimental values for HI data points; we see on 
this figure the good agreement of the QCD fit with the data points, which supports the validity of the description of 
the Pomeron in terms of partons following QCD dynamics. In Fig. 3 we also compare the scaling violations obtained 
in the fit to the HI measurement. 



5 If one fits the values of ap and %, one gets the same values as in Section II. 

6 Fit 2 shows a higher \ 2 compared to fit 1 (the \ 2 P er degree of freedom is 1.15 for fit 1 and 1.25 for fit 2). Furthermore, fit 
2 is quite unstable: changing the parameters a little modifies the gluon distribution at high z, and removing some data points 
at low or medium (3 for instance destroys this fit completely. The quark distribution is similar for both fits, but the gluon 
distribution tends to be quite different at high values of z. This can be easily explained as no data above z — 0.65 are included 
in the fits. Thus there is no constraint from the data at high z. Therefore, we keep fit 1 as the best fit for HI diffractive 
cross-sections. 
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Let us discuss the assumptions we made in the QCD fits. We took Q 2 , = 3 GeV 2 and we noticed that the x 2 of 
the fit increases significantly if we try to lower this value of Qq (x 2 /dof=1.5 if <3o=l GeV 2 ), indicating that the fit 
cannot be extended to lower Q 2 . We will come back to this point in the following. 

We performed the fits including three Chebyshev polynomials. Including one more polynomial does not improve 
the x 2 of the fit, while two polynomials are also not enough to get a good convergence of the fit (x 2 /dof=1.7). 

The paramctrizations used in this analysis are given in formulae (3) and (4) : indeed, the behaviour at large z is 
forced to go to when z goes to 1 by the term exp(-^j) which includes an essential singularity at z = 1. We tried to 
perform the analysis with a = imposing the behaviour at large z with a term in (1 — z) a . With these assumptions, 
we did not get any fits with a good x 2 (x 2 /d° r >2). The form of the parametrization for the parton distribution at 
large z plays a non-negligible role in the fit procedure to ensure a quick convergence at high z. However, this term is 
large in a region where no data are included in the fit. On a more physical ground, it is interesting to note that this 
behaviour is in contrast with the one for proton (or even pion, cf. Drell-Yan analysis) structure functions for which 
the behaviour of the distribution tail is seen in a rather sizable region of z. 

The positive scaling violations of diffractive structure functions lead to an important gluon content zG[z) in the 
Pomeron (even at large z), which is what we get from QCD fits (see Fig. 3). 

Moreover, with this gluon distribution, larger scaling violations at low (3 are naturally associated with larger values 
of gluon splitting functions when [3 gets lower. Now, it is also clear that, due to large measurement uncertainties, 
scaling violations and thus the determination of a s zG are not very precise. The uncertainty on the value of a s is 
negligible in this analysis, and we can get direct access to the uncertainty on the gluon density (determined from 
QCD fits). We get from the investigation of the scaling violations 

^25%. 

We believe that this (even not yet precise) error determination is a relevant piece of information, showing that the 
gluon density is not yet very well constrained. 

In conclusion, the gluon content of the Pomeron is confirmed to be large (w.r.t. [1]). The flat gluon solution 
is favoured, and another result (fit 2) appears if we allow that the parton distributions can possibly vanish in the 
kinematical range. Moreover this second solution is disfavoured due to its worse x 2 an d instability and thus, we 
will not use this second solution in the following. The accuracy of the gluon determination is still quite poor. This 
provides a very nice perspective for diffractive analyses of more recent data where more accurate results will lead to 
a more precise QCD analysis. 
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FIG. 1. Quark flavour singlet (zS, left) and gluon (zG, right) distributions of the Pomeron deduced as a function of z, the 
fractional momentum of the Pomeron carried by the struck parton, from the fit on HI data points with Q 2 > 3 GeV 2 . The solid 
(dashed) curve shows fit 1 (fit 2) as discussed in the text (x 2 /dof = 177.1/154 = 1.15 for fit 1 and x 2 /dof = 192.5/154 = 1.25 
for fit 2 with statistical errors only). The parton densities are normalised to represent xp times the true parton densities 
multiplied by the flux factor at xp = 0.005 (this will always be the case in all the following figures). 
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FIG. 2. The HI data points for x F F° (3) are shown with the result of the QCD fit (fit 1) described in the text; the result of 
the fit is shown only in bins included in the minimization procedure. 
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FIG. 3. The structure function x-p ■ at xp = 0.005, presented as a function of Q 2 in bins of (3, over the full Q 2 range 

accessed with the HI data. The superimposed line represents the result of the QCD fit to HI data described in section IIIA. 



B. Parton distributions from ZEUS data 



We have redone this QCD analysis with ZEUS diffractive cross-section measurements [2] applying the same cuts. 
Note that for the ZEUS fits, there are only 6 parameters since the method of measurement is not sensitive to the 
Reggeon contribution [2]. The Pomeron intercept is fixed to ap = 1.13 as determined from Regge fits on ZEUS data 
[2]- 

The resulting parton distributions are presented in Fig. 4 together with the HI results. As the Pomeron intercepts 
are different for HI and ZEUS data, the parton densities are normalized to represent xp times the true parton densities 
multiplied by the flux factor at xp = 0.005. 
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We notice some large discrepancies, especially in the gluon density 7 . The gluon density derived from ZEUS 
measurements is lower that the HI fit results. The data are compared directly in Fig. 6 where we rescaled the 
HI data to the ZEUS bins (the rescaling is not very model-dependent since the shifts in Q 2 and (3 between both 
experiments are quite small, and we performed the rescaling using models in Ref. [16] [17], which lead to the same 
result). In addition to this treatment of data, about 10% (30 %) of the cross-section should be subtracted to take 
into account proton dissociation contribution in the HI (ZEUS) cross-section measurement. From Fig. 6, we notice 
that the main differences are located in the highest Q 2 bins (Q 2 ~ 60. GeV 2 ) after this additional correction and 
could be due to the different methods used to do this measurement (rapidity gap method or so-called Mx method 
[1,2]). However, there are only 30 ZEUS data points entering in the QCD analysis compared to 161 for HI analysis, 
the experimental uncertainties being of the same order. Thus, the statistical power of the ZEUS fits is considerably 
lower. We can estimate the statistical uncertainty on the gluon distribution to be about 50% for ZEUS data (to be 
compared with 25% for HI data). The differences in the quark densities are much smaller taking into account the 
fact that there is no ZEUS data at Q 2 = 4.5 GeV 2 and low (3. Differences in the parton distributions are clearly due 
to small differences in the data (see Fig. 4). We see differences in the Q 2 =60 GeV 2 , 0=0.7, or Q 2 =8 GeV 2 , (3=0.2 
bins. 

As the fit results are quite different for both experiments, we choose to make predictions for the longitudinal and 
charm structure functions for both sets of parameters. 



7 Note that such different gluon densities can lead to similar values of the structure functions (see Fig. 4) because the 

diffractive structure function is more sensitive to the quark density (the photon couples directly to quark). This will however 
lead to differences in the charm structure function as it is much more sensitive to the gluon structure function. 
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FIG. 4. Quark flavour singlet (zS, left) and gluon (zG, right) distributions of the Pomeron derived from HI diffractive data 
(x 2 /dof = 177.1/154 = 1.15, full lines) and ZEUS diffractive data {x 2 /dof = 29.3/(30 - 6) = 1.22, dashed lines). 
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FIG. 5. The ZEUS data points on x w F° (3) are shown with the result of the QCD fit described in the text. 
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FIG. 6. Comparison of HI and ZEUS data on xpF®^ (respectively black squares and open triangles) with the result of the 
QCD fits: full line for the HI QCD fit result and dashed line for ZEUS. Note that the HI points have been rebinned to ZEUS 
Q 2 and f3 values. 
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IV. GLOBAL FITS 



A. Higher-twists and global fits 



The HI fit (fit 1) is our reference and we attempt to refine it by extending the number of data points included 
into the fit procedure. The problem is that the cuts described above are essential to avoid the experimental domain 
where non-perturbative power contributions could be large and we cannot suppress these cuts without subtracting 
in another way these power contributions (or higher- twist terms). However, some parametrizations for higher- twist 
contributions have been proposed [16,17] after a dedicated analysis of HI data. 

We use the results of Ref. [16,17] to subtract the higher twist component from the measurements at high (3 

and we redo the QCD analysis with all data points included, which gives 224 data points for HI. 

The resulting parton distribution are presented in Fig. 7 for the fit to the HI data using the higher twist deter- 
mination of Ref. [16] 8 . We notice the perfect agreement between parton distributions extracted with this procedure 
and with our former analysis described above. However, the statistical power of the new procedure is better due 
to the larger number of degrees of freedom. We also note that the x 2 /dof is smaller when we perform the fit after 
subtracting the higher twist component (1.05 compared to 1.15, see Fig. 7). This implies that a small residual higher 
twist component is still present near /3=0.65 which influences the quality of the fit, and thus, it is quite dangerous to 
cut directly in the data. 

As we now include all data in the fit, the assumption we made on the a parameter of Equation (4) may play a 
role. We checked as before that the convergence is not possible if we include a polynomial instead of an exponential 
form (the \ 2 P er degree of freedom increases up to 2.). However, the value of a is not fundamental as another value 
(a = 10~ 3 ) also leads to a good fit. 

The parton distributions extracted with the global procedure are now our reference distributions. The parameters 
of the gluon and quark densities are given in appendix A2. 

From the ZEUS data, we also get comparable results with or without higher twist corrections, leading to the same 
parton distributions as shown in Fig. 4. 
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FIG. 7. Parton distributions extracted from HI data with the standard QCD fits (\ 2 /dof — 177.1/154 = 1.15, full lines, fit 
A) compared to distributions extracted from the global fit analysis using all HI measurements (x 2 /dof = 228.4/(224-6) = 1.05, 
dashed lines, fit B). 



8 We checked that using another higher twist parametrisation from Ref. [17] which uses a BFKL-dipole type approach [18,19] 
does not change our results. 
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B. Comparison of the gluon density in the proton and in the Pomeron 



As can be seen in Fig. 8, the gluon density derived from our global QCD fits of HI diffractive measurements is very 
different from the one derived from QCD fits of F2 [24]. However, since the uncertainty of our determination is large 
(77- > 25% as we have mentioned above), one can wonder whether could be described using the same gluon 

density as for the proton total structure function F2. 

We have taken the gluon density of Ref. [24] (we simply identify Bjorken-rr with (3 from this parametrisation) and 

we have redone the QCD analysis of F®^ HI data by fitting only the sea quark distribution using directly the 
gluon density inside the proton. As explained above we have also used the prescription of Ref. [16] to subtract the 
higher-twist contribution from diffractive measurements, which gives 224 data points and 3 parameters in the fit. 
This procedure leads to "X 2 /dof = 304.2/(224 — 3) = 1.37 which is worse than our reference fit for diffractive data 
(x 2 /dof — 228.4/(224 — 6) = 1.05). This result suggests that the gluon density of the Pomeron is quite different from 
the gluon density of the proton 9 . The difference is smaller between the result of the ZEUS fit and the gluon density 
inside the proton (in the case of the ZEUS fit, the gluon density in the Pomeron is found to be smaller than in the 
HI case). Both HI and ZEUS gluon densities in the Pomeron are however incompatible with the gluon in the proton 
at high x or fj. 

The fit results compared to the data show large discrepencies at low z in all Q 2 bins. It is impossible to describe 
diffractive data at low (3 using the same gluon density as in the proton. 
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FIG. 8. Gluon density fitted from diffractive data (global fits, full lines) compared with the gluon density extracted from F2 
QCD fits (dashed lines). 



9 Note that we have also redone this procedure with a free normalization for the gluon distribution. However, this normalization 
is not changed during the fitting procedure. Thus we get the same results as before. 
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C. Fit extrapolations in Q 2 



We now extrapolate our QCD fits using HI and ZEUS data to low and high Q 2 to look for a possible kinematical 
domain where the differences in the HI and ZEUS parton distributions can lead to differences in . 

As we mentioned before in section IIIA, it is difficult to lower the starting scale Q\ to get a prediction for 
at low Q 2 because the quality of the fit is much worse. Thus, we perform a low Q 2 backward evolution using Qq=3 
GeV 2 . On the other hand, it has been recently noted by the HI and ZEUS collaborations [20] that the HI 1995 and 
ZEUS 1995 preliminary data at low Q 2 cannot be described using our perturbative QCD approaches, which might 
be an indication that perturbative QCD or a DGLAP fit to F 2 D is not valid in this region of phase space without 
saturation effects [21]. 

At high Q 2 (see Fig. 9) the differences between both extrapolations remain small as was the case at medium Q 2 , 
except at high xp, which is an artefact of the absence of Reggeon contribution in the ZEUS fit. 

The fact that no large difference is found between both extrapolations whereas we have large differences in the 
gluon density shows that we need other measurements more sensitive to the gluon density, like the charm component 
to F 2 D , or some direct final state studies (jet production, multiplicity studies [23]). 
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D. Charm contribution and longitudinal diffractive structure functions 



From global QCD fits of the diffractive structure function F 2 , we have derived reference parton distributions 
(see Fig. 7). With these distributions, we can now give some predictions for the longitudinal and charm structure 
functions in diffraction. 

We present values for the longitudinal structure function in Fig. 10, which is the sum of two components, namely 
the perturbative part coming from the QCD fit and the higher-twist component. The higher-twist component of the 
longitudinal structure function becomes dominant at high value of f3 whereas the low (3 domain is dominated by the 
perturbative part of Fl deduced from QCD fits. We can also evaluate values of R = ^ in diffraction, which is 
illustrated in Fig. 11. Here again, the high f3 domain represents the higher- twist component. 

We note that the differences between the HI and ZEUS fit results are quite small in F^. However, a measurement 
of F£ or R would be quite interesting to constrain the higher twist component at high f3. The R values are more 
different between HI and ZEUS predictions, but the error bars on these extrapolations are quite large. 

The charm predictions are shown in Fig. 12, 13. We note the large differences between the ZEUS and HI fit results 
which shows the importance of this direct measurement, which is directly sensitive to the gluon density. The charm 
contribution to F® is expected to be larger than for the proton structure function (about 40% at (3 ~0.1, Q 2 ~ 10 
GeV 2 ) if we consider the results of the HI fit, whereas the results of the ZEUS fit are more similar to the proton 
values 10 . 




FIG. 10. Prediction for the longitudinal structure function xp ■ (at fixed xp = 0.005) presented as a function of z = (3 

for different values of Q 2 (HI: full line, ZEUS: dashed line). 



10 Preliminary HI and ZEUS measurements [22] show that the HI (resp. ZEUS) charm cross-section is closer to the results of 
the predictions coming from the ZEUS (resp. HI) QCD fit. More precise data are needed to solve these discrepencies. 
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V. EXTRAPOLATION TO TEVATRON AND COMPARISON WITH CDF DATA 



The QCD fits we obtained from HERA data allow us to make direct comparisons for measurements at the Tevatron. 
It is quite interesting to be able to directly test factorization breaking between HERA and the Tevatron using the 
measurements performed at both accelerators. We thus compare the extrapolations of the HI and ZEUS QCD fits 
to the recent CDF diffractive jet cross-section measurement [7]. The result is given in Fig. 14. We note a large 
discrepancy both in shape and normalization between HI predictions and CDF data, clearly showing factorization 
breaking. However, the ZEUS fits are more compatible in normalization with the CDF measurement even if the shape 
is not described properly. We know that the ZEUS gluon density is between 2 and 3 times smaller than the HI one. 
The predictions for the Tevatron are thus expected to be very much different by a factor 2 to 3 between HI and ZEUS 
since they correspond to single Pomeron exchange. If the large statistical and systematic uncertainties on the gluon 
density are taken into account (about 50% for ZEUS, 25% for HI), ZEUS data are compatible with factorization at 
low (3. One should, however, question whether one has the right to extrapolate ZEUS results without introducing an 
additional Reggeon component. Namely the CDF measurement is in a region in xp where the Reggeon contribution is 
important, in contrary to the ZEUS measurement. The error bar on the ZEUS fit extrapolation is thus enhanced due 
to this uncertainty. Concerning the extrapolation of the HI fit and the comparison with CDF data, one must notice 
that the CDF data lie primilarly at low values of /3 where there are few HI data points. Moreover, in this kinematical 
domain, the HI data have the tendency to lie more in the high xp region where the Reggeon contribution is important 
and not well constrained by the fit. Thus, the extrapolation to the CDF domain suffers from large uncertainties. A 
combined fit using CDF data to constrain the low [5 region and the HERA data to constrain the high (3 domain would 
thus be of great interest. 

More precise data from HERA and more detailed comparisons with Tevatron are thus needed to study precisely 
factorization breaking between both experiments. The discussion of eventual higher twist contributions is clearly also 
valuable in order to reach conclusions on the shape at larger /?. It is thus important to get an accurate measurement 
of the gluon density in the Pomeron at low values of (3 from HERA. Furthermore the Forward Proton Detector [25] 
installed by the DO collaboration for Run II will be of great help to get a direct measurement of the diffractive 
structure functions. 
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FIG. 14. CDF data j3 distributions compared with extrapolations from the parton densities of the proton extracted from 
diffractive measurements from HI and ZEUS collaborations 
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VI. CONCLUSION 



Let us again go through the points discussed in the introduction: 

-We have developped a new method to get the Pomeron and Reggeon exponents using Regge fits, confirming a 
higher value of the Pomeron intercept in "hard" diffraction than in soft processes using a more precise method. 

-We propose (and check the stability) of a new NLO DGLAP fit to the diffractive structure functions. We have 
obtained parton distributions for the Pomeron using HI and ZEUS data. 

-HI and ZEUS data both require a large gluon component inside the Pomeron but the gluon densities are quite 
different. The quark densities are found to be similar in shape and magnitude. Using a paramctrisation for the 
higher twist contribution to allows us to get a description of the full set of data using the NLO DGLAP evolu- 
tion equations, and leads to the same parton distribution as before. The longitudinal diffractive structure function 
measurement could allow us to further constrain these higher twist paramctrisations. 

-Differences in the Hf and ZEUS gluon distributions lead to significant differences in the longitudinal and charm 
contributions to F® . The low Q 2 extrapolation shows some scaling violations which do not seem to be present in the 
data, which might be a sign of saturation effects. 

-The same difference leads to diverging extrapolations for single diffraction at the Tevatron. While the Hl-based 
extrapolation leads to a strong factorization breaking in both shape and normalization (at least a factor of 10 with 
respect to CDF data), the ZEUS-based extrapolation leads to a possible compatibility with factorization at least in 
the small [3 (large diffractive mass) range. 

All these questions deserve a development of experimental and theoretical work in this domain, with the prospects 
of undestanding key aspects of QCD properties of "hard" diffraction. The possibility of diffractive production of 
heavy objects (jets, Higgses) at tevatron even enhances the interest of such studies. 
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Appendix Al: Decorrelation method to get the Pomeron and Reggeon intercepts 



Rewriting equation (1) as 

F 2 D(3 \Qlp k ,x P:i ) = fp/pixp.OA^QlPk) + h </v (^ :i )B jk (Qlf3 k ) (5) 

A fit is then performed in each (<3|,/3fe) bins. A and B are only related to the Q?,(3 k bin considered, which means 
that they can be fitted within this bin by minimizing the following XjkiQji Pk) functions 
We have X 2 = T,j E fc Xjk(Q 2 j,Pk) and 

2 fn2 R sr {h/ v {^MAQ 2 rPk)+fu/ v {^s)B {Q],f3 k ) - F 2 D{3) ' exper -(Qlp k ,x FA ))i 
X3kiQ ^ k) \ A ^(3).^r. {Q% p M)Y (6) 

where the i index runs over the number of points with different values of xp in the Q?,f3 k bin considered and 
^pD(3),exper. . g statistical error on F 2 D (denoted a-i in the following). We use the notations 

fp/p{xp,i) = /p,i (7) 
fu/ P ( x R,i) = fn,i (8) 
F 2 D ^ exper -(Qlp k ,x F , i ) = F i . (9) 

As Aj k and Bj k do not depend on ipj, they can be calculated directly by solving the following system 

dx 2 (QlPk) 



dA ]k 
d X 2 (Q 2 ,[3k) 



dB 



= (10) 
= (11) 



which gives 



7T" 



^ k ~ * 2f 2 *2f 2 '„ , * (^) 

E ^/e iFj y-v ^/jp t 2/p jfi ^/R,i/p,i 

Thus, for each Q^,(3 k bin, A,-fc and Bj k are well defined values and depend on ap(0) and «r(0) via the flux factors. 

In conclusion, all Aj k and Bj k parameters should not be introduced without care as free parameters in the global 
Regge fit as they satisfy the two last expressions in terms of Fi, /p t ; and /r,;, which reduces considerably the number 
of free parameters. Indeed only 2 free parameters remain : cnp(0) and «r(0). 
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Appendix A2: ZEUS and HI fit parameters 



can be expressed as a sum of two factorized contributions corresponding to a Pomeron and secondary Reggeon 
trajectories. The Pomeron and Reggeon fluxes are assumed to follow a Regge behaviour with linear trajectories 
oep,u(t) = oep,u(Q) + "p.r^j such that 



fw/p,R/p( X TP) 



2«ip,ir ,(t)-l 



dt 



(14) 



x 



p 



where \t min \ is the minimum kinematically allowed value of |t| and t cut = — 1 GeV is the limit of the measurement. 



? D(3) 



(Q , P, xp) = f P/p (x F )F^(Q 2 , p) + N n f R/p (x P )F^(Q 2 , (3) . 



The values of Q;p j r(0) are free parameters while -Brr and a P R are taken from hadron-hadron data (a p = 0.26 

GeV -2 , a n = 0.90 GeV -2 , B F = 4.6 GeV -2 , B F = 2.0 GeV -2 The Pomeron intercept coming from Regge fits is 
a P (0) = 1.20 ± 0.02 for HI, and a P (0) = 1.13 ± 0.04 for ZEUS. 

The Reggeon structure function is needed only to describe HI data, and the value of Qr(0) = 0.62 ± 0.02, with 
the normalization N u = 14.2 ±0.1. The values of F 2 R (g 2 , (3) are taken from a parameterisation of the pion structure 
function [9]. 

The parametrisations for the gluon and quark densities (zS q (z, Q 2 ) = u + u + d + d + s + s) in the Pomeron are 
given at Ql=3 GeV 2 : 



zS(z,Q 2 = Q 2 ) = 
zG(z,Q 2 = Ql) = 



Y^cf-P^z-l) 



e *-i 



J2q G) -p 3 (2z-i) 

3=1 



(15) 
(16) 



where Pj{0 is the j member in a set of Chebyshev polynomials, which are chosen such that Pi = 1, P 2 = C an d 

P j+1 (C) = 2CP J -(C)-P j -i(C)- 



parameters 


HI 


ZEUS 




0.18 ± 0.05 


0.41 ±0.02 




0.07 ± 0.02 


-0.16 ± 0.03 




-0.13 ± 0.02 


-0.11 ± 0.02 




0.82 ± 0.40 


0.53 ± 0.30 




0.22 ± 0.06 


0.28 ± 0.25 


cP 


0.01 ± 0.04 


0.02 ± 0.11 



Table I- Pomeron quark and gluon densities parameters 
The parametrisation for the higher-twist contribution to the diffractive structure function is the following 



7 D(3),III 



( [ .Li I ( Qo.\ ( i~f Q 2 



xp 



Q 2 



ln( 



AQ 2 (3 



where 



"2,4 = nl + n\ln 
with the following values for the parameters 



, ,Q 2 



1.75) /3 3 (l-2/3) 2 



(17) 



(18) 
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parameters 


HI 


ZEUS 


C 


0.035 ± 0.016 


357. ± 106. 


x 


0.40 ± 0.02 


0.0008 ± 0.0002 


nl 


1.43 ± 0.08 


1.00 ± 0.10 


n\ 


0.00 ± 0.05 


0.37 ± 0.10 



Table II: Parameters obtained for the higher-twist parametrisation 
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